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Abstract 

This paper deals with the derivation and analysis of the the Hall Magneto- 
Hydrodynamic equations. We first provide a derivation of this system from a two- 
fluids Euler-Maxwell system for electrons and ions, through a set of scaling limits. 
We also propose a kinetic formulation for the Hall-MHD equations which contains as 
fluid closure different variants of the Hall-MHD model. Then, we prove the existence 
of global weak solutions for the incompressible viscous resistive Hall-MHD model. 
We use the particular structure of the Hall term which has zero contribution to the 
energy identity. Finally, we discuss particular solutions in the form of axisymmetric 
purely swirling magnetic fields and propose some regularization of the Hall equation. 
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1 Introduction 

This paper deals with the derivation and analysis of the the Hall Magneto- Hydrodynamic 
(Hall-MHD) equations. The Hall-MHD model currently receives an increasing attention 
from plasma physicists. It is believed to be the key for understanding the problem of 
magnetic reconnection. Indeed, space plasma observations provide strong evidence for the 
existence of frequent and fast changes in the topology of magnetic field lines, associated 
to violent events such as solar flares [11] . However, magnetic reconnection cannot be 
described in the framework of ideal MHD, due to the frozen-field effect. Indeed, in ideal 
MHD, due to the Faraday equation and ideal Ohm's law, the magnetic field is essentially 
passively transported by the fluid velocity. Therefore, the topology of the magnetic field 
is preserved, even in the magnetic field lines are deformed by the flow. In order to 
break this passive magnetic field transport by the fluid flow, one is led to re-introduce 
the Hall terms which was neglected in ideal MHD. In spite of its increasing importance 
for physical applications, the Hall-MHD model has received very little attention from the 
theoretical viewpoint (see e.g. [13 Ej) and the purpose of this paper is mainly to propose 
a framework for the derivation and analysis of the Hall-MHD problem. 

We first provide a derivation of this system from a two-fluids isothermal Euler-Maxwell 
system for electrons and ions, through a set of scaling limits. The two-fluids model of 
plasma is known as Braginskii model [2] and was justified in [9] (see also [5]) on the 
basis of a concurrent hydrodynamic and zero electron to ion mass ratio limit. As usual 
in MHD models, zero electron to ion mass ratio limit, zero Debye length limit and zero 
displacement current limits have to be taken. Then, the main point is to examine the 
orders of magnitude of the various terms arising in the generalized Ohm law (which is the 
electron momentum equation with zero inertia) and in the current equation. Basically, 
Hall MHD is obtained when the electron and ion velocities have difference of order unity 
and when this difference is introduced inside the generalized Ohm law. 

Then, we propose a kinetic formulation for the Hall-MHD equations which contains as 
fluid closure different variants of the Hall-MHD model. The kinetic formulation consists 
of a Fokker-Planck equation for the ions and a set of fluid equations for the electrons cou- 
pled through quasineutrality. The Fokker-Planck operator models electron-ion collisions 
and contributes to relaxation of the velocities and the temperatures of both species to a 
common value. This kinetic model was in particular justified in [TJ. The MHD equations 
are obtained by taking the fluid moments of the ion Fokker-Planck equation and closing 
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the resulting equations by a Maxwellian assumption. The resulting two-temperature re- 
sistive Hall MHD model consists of conservation equations for the density, momentum, 
energy and magnetic field combined with an evolution equation for the electron temper- 
ature and with the generalized Ohm law. The proposed hybrid ion-kinetic, electron-fluid 
model bear strong analogies with models used in the literature for numerical simulations 
such as [22], but the purpose is to highlight its mathematical structures. Indeed, a par- 
ticularly interesting special case is when the electron and ion temperatures are equal to 
the same constant value (isothermal single temperature resistive Hall MHD). In this case, 
we can rephrase its kinetic formulation in the form of a coupled Fokker-Planck Faraday 
system, which exhibits an entropy dissipation identity. Surprisingly enough, the kinetic 
formulation of standard ideal MHD is deduced by neglecting the Hall term in Faraday's 
equation but keeping it in the kinetic equation. 

The theoretical analysis focuses on the existence of global weak solutions for the in- 
compressible viscous resistive Hall-MHD model written as follows: 



d t u + u-Vu + Vp = (V x B) x B + Au, (1.1) 

V-u = 0, (1.2) 

d t B - V x (u x B) + V x ((V x B) x B) = AB, (1.3) 

V-B = 0, (1.4) 



where u(x,t) and B(x,t) are the fluid velocity and magnetic field, depending on the 
spatial position x and the time t. The result is valid on a square domain Q of M 3 with 
periodic boundary conditions. We stress that the important contribution of this work is 
the account of the last term of the left-hand side of (ll.3|) . known as the Hall effect term. 
The main theorem of this work is stated as follows: 

Theorem 1.1 LetVt = [0, l] 3 . Assume thatu G (L 2 (fi)) 3 ; B G (L 2 (tt)) 3 withV '-u = 0, 
V ■ Bq = 0. Then, there exists a global weak solution (u,B) for the Hall MHD problem 
/ TO) , {pp. Moreover, we have (u,B) G L°°((0, T), L 2 (fi)) n L 2 ((0,T), H\n)) and 
d t u G Ls((0,T), J ff- 1 (fi)) ; d t B G Lt((0,T),#- 2 (fi)). Additionally, the following energy 
inequality holds: 

j t S(t) + \\Vu\\l 2{a) + \\VB\\ 2 L2{Q) < 0, (1.5) 

with 

£(t) = l\\u\\ 2 L2{n) + \\B\\ 2 L2{n) . (1.6) 



In fact, the main difficulty of this work is concentrated in the treatment of the Hall 
term. So, we will show a preliminary result for the following Hall problem 

d t B + V x ((V x B) x B) = AB, (1.7) 
V-5 = 0, (1.8) 
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and will provide a detailed proof. The proof uses the particular structure of the Hall term 
which has zero contribution to the energy identity. The proof of the existence for the 
coupled system is then a direct consequence of the energy inequality (11.51) and will only 
be sketched. This result is up to our knowledge the first theoretical result for Hall MHD. 
Global existence for standard viscous resistive incompressible MHD has been previously 
proved by Duvaut & Lions [9]. A bifurcation analysis of the Hall-MHD problem in view 
of the question of magnetic reconnection is performed in [T2]. Numerical methods for 
solving the Hall-MHD problems can be found e.g. in [IJ HJ EH EES] 

Remark 1.1 In a general domain Q, the physically relevant boundary condition is the 
perfectly conducting wall boundary condition, which consists in assuming zero normal 
component of the B field and zero tangential component of the electric field. Here, the 
electric field is the quantity inside the curl operator, namely 



[ n x (V x B) + n x ((V x B) x B) = on dVt. 
From these conditions, we deduce that 

B ■ (V x B) = 0, 

which means that there is no helicity on the boundary. Because of the nonlinearity of this 
boundary condition, the methods developed below do not apply. 

Finally, we discuss particular solutions of the Hall problem in the form of axisymmet- 
ric purely swirling magnetic fields and propose some regularization of the Hall equation. 
For axisymmetric purely swirling magnetic fields, the Hall problem reduces to a viscous 
Burger's equation. By neglecting the resistivity, the resulting inviscid Burger's equation 
shows shock wave solutions which are known in physics textbooks as KMC waves for 
Kingsep, Mokhov and Chukbar [IB] (see also [21 E])- They only exist if the Hall term 
is present. This effect also generates boundary layers which lead to nonlinear boundary 
conditions (see e.g. [20]). Focusing on the non-resistive Hall problem itself, we propose 
a regularization consisting in restoring the displacement current in the Ampere equa- 
tion. We then provide two equivalent formulations of this regularized problem which are 
obtained when either the current or the electric fields are eliminated from the system. 

The organization of the paper is as follows. In Section [21 we propose a derivation 
of the Hall-MHD model from the two-fluids Euler-Maxwell model under suitable scaling 
hypotheses. Then, Section [3] is devoted to the presentation of the kinetic formulation 
of the Hall-MHD problem. Section H] is focuses on the proof of the existence of global 
weak solutions for the incompressible viscous resistive Hall-MHD equations. Section [5] 
discusses the particular case of axisymmetric purely swirling magnetic fields and proposes 
a regularization of the Hall problem by means of a re-introduction of the displacement 
current in the Ampere equation. Finally, a conclusion is drawn in Section [6j 



E = (V x B) x B + V x B. 
This leads to the nonlinear boundary conditions 




on dVl, 
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2 Derivation of the Hall-MHD equations from fluid 
equations 

In this section, we briefly motivate the derivation of the model we are considering. For 
simplicity, we consider the compressible inviscid model and later on change to viscous 
incompressible flow. We start from the two-fluid isothermal Euler-Maxwell system for the 
electrons and ions, where we assume that the electron and ion temperatures are the same 
given constant: 



c~ 2 d t E - V x £? = -/i j, 

e V ■ E = p, 

d t B + V x E = , 

V-B = 0, 

p = e(rii - n e ), 

j = e(riiUi - n e u e ). 



where n e and rij are the electron and ion densities, u e and Ui, the velocities, T, the common 
electron and ion temperature, m e and the masses, e denotes the elementary positive 
charge, and we assume singly charged positive ions, rj is the resistivity due to the electron- 
ion collisions. E, B, p, j are respectively the electric field, the magnetic field, the charge 
density and the current density. e , Po and c are respectively the vacuum permittivity, 
the vacuum permeability and the speed of light, related by the relation e^p^c 2 = 1. For 
simplicity, we assume monoatomic perfect gas equations of states for both the electrons 
and ions. We make the Boltzmann constant equal to unity which means that we measure 
temperatures in units of energy. The last terms at the right-hand sides of the second 
and fourth equations are the contributions of the electron-ion collisions to the momentum 
equation of each species. The two terms sum up to zero which expresses the conservation 
of total momentum in such collisions. 

We introduce scaling units n , uq, E , B , Xq, to, po, jo for respectively the densities, 
velocities, electric field, magnetic field, space, time, charge and current. We assume that 
these units are related by the following relations: 



d t n e + V ■ {n e u e ) = 0, 



(2.1) 



m e (d t (n e u e ) + V(n e u e <g> u e )) + V(n e T) = 

= —en e (E + u e x B) — e 2 rjn e rii{u e — itj), 



(2.2) 



d t rii + V ■ {riiUi) = 0, 

mi(d t {niUi) + WiriiUi <g> Ui)) + V(njT) = 



= erii(E + Ui x B) — e 2 r\nin e {ui — u e ), 
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The first relation means that we observe the system at the convection time scale. The 
second relation states that the convection velocity is that of the ion thermal speed. The 
third relation is typical of a MHD scaling and states that the main contribution to the 
electric field is induction due to the motion of the charged fluid. Finally, the last relation 
expresses the consistency between the density and charge units. 
Then, six dimensionless parameters appear: 

2 m e 2 eE x e 2 r]n u x 
mi 1 1 

7=—, A =— 2' V = • 

which have the following interpretation, e 2 is the electron to ion mass ratio and is very 
small, a 2 is the ratio of the electric energy to the thermal energy. (3 measures the 
relaxation frequency of the electron and ion velocities due to collisions. 7 is the ratio of 
the fluid velocity to the speed of light. A is the scaled Debye length and measures the 
closedness to quasi-neutrality. 77 is the ratio of the charge current scale to the electron or 
ion current scales. Since the charge current is the difference of these two particle currents, 
it may be much smaller than any of them due to charge neutrality. Therefore, the scale 
ratio 77 may be either 0(1) or <C 1 according to the situations. 
The dimensionless two-fluids Euler-Maxwell system is written: 

d t n e + V • {n e u e ) = 0, 

e 2 (d t (n e u e ) + V{n e u e ® u e )) + V(n e T) = 

= —a 2 n e (E + u e x B) — /3n e rii(u e — ui), 

d t rii + V • (riiiii) = 0, 

d t {n,iUi) + V(niUi <g> it*) + V{n,{T) = 

= a 2 rii(E + Ui x B) - f3rnn e {ui - u e ), 
2 

7 2 d t E - V x B = -^U 



ft" 



! A 2 ^ 



ft 2 A 2 V -E = p, 
d t B + Vx£ = 0, 
V-S = 0, 
p — rii — n e , 

j = -(riiUi - n e u e ). 
V 

The compressible MHD equation corresponds to the simultaneous independent four 
limits 

1. e 2 — > 0: this corresponds to the neglect of the convection term in the electron 
momentum equation. The resulting equation is usually referred to as the generalized 
Ohm's law. 
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2. A 2 — > 0. This gives rise to quasineutrality, i.e. the fact that the local electron and 
ion densities are everywhere the same. We now denote by n their common value: 
n e = rii = n. 

3. 7 2 — > while keeping = 1. This leads to the neglect of the displacement current 
in Ampere's equation and gives rise to the standard magnetostatic Ampere law. 

The resulting system is the so-called compressible isothermal resistive Hall-MHD equa- 
tions. Denoting by u the ion velocity, this system is written 

d t n + V • (nu) = 0, 

d t {nu) + V{nu ®u) + V(2nT) = a 2 7] j x B, (2.3) 

VxB = j, 

d t B + V x E = 0, 

V-5 = 0, 

j = -n(u-u e ), (2.4) 
V 

T 7 x B 3n 

E + uxB = -V(ln n) + r) J - + ^j, (2.5) 

where we highlight the momentum conservation eq. (12.31) . the current equation (12.41) and 
generalized Ohm's law (12. 5p . Note that the T Vlnn term at the right-hand side of (12. 5p 
has no contribution since the curl operator in the Faraday equation cancels it. However, 
this cancellation is no more true in the general gas dynamics case because may not 
be a gradient in general, where p e is the electron pressure. 

In all what follows, we assume a 2 rj = 1 in order to keep the Lorentz force term in (12.31) 
of order 1 in all the various scalings below. Then, eqs (12.31) . (12 .4p and (12. 5p are written : 

d t {nu) + V{nu ® u) + V{2nT) =jxB, (2.6) 
\j = n{u-u e ), (2.7) 



a 

E + U X B = At 

a 2 



-TV (Inn) -' 



+ —J, {2.1 

a 4 



the other equations being unchanged. There are only two dimensionless parameters left: 
^2 and ^ and they only appear in (12 .7p and in (12.81) . So, the various types of MHD model 
correspond to the various choices of scalings for these two parameters. In particular, we 
have 

1. If both \ — y and -4- — > 0, then the generalized Ohm's law reduces to the standard 
ideal Ohm's law while the electron and ion velocities become identical: 

E + u x B = 0, u e = u. 

This yields ideal MHD. 
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4 ->• but 4 
but the electron and ion velocities are still identical: 



If 4 — ► but — >• 1, then the resistive term in the generalized Ohm's law is kept 



E + u x B = j, u e = u. 
This gives rise to resistive MHD. 

3. If 4 -> 1 but 4 -»■ 0, then, the ion and electron velocities differ and additionally, 
the generalized Ohm's law has the form: 

7 X B 

E + ux B = -TV(lnn) + . 

n 

As already mentioned, the first term at the right-hand side has no contribution. 
The second one is the Hall term. This gives rise to the Hall MHD. 

4. Finally, if both 4 — > 1, 4 — >• 1, then, the ion and electron velocities differ and both 
the Hall and resistive terms appear. 

j x B 

E + ux B = -TV (Inn) + + j. 

n 

Our study takes place in the context of the last regime, where both the resistive and 
Hall terms are equally important. Additionaly, we assume incompressible viscous fluid 
motion. In this case, the Hall MHD system can be written according to (ll.ll) - (11.4j) . 

We note that it is easy to extend this system to the viscous isentropic resistive com- 
pressible Hall MHD as follows (assuming all the physical constants equal to 1 except the 
viscosity here denoted by v): 

d t n + V • (nu) = 0, (2.9) 
d t (nu) + V(mt <g> u) + Vp(n) = (V x B) x B + z/V • (Vit + (Vu) T ), (2.10) 

d t B + V x (b x u + ^ V X B ^ X B \ = -V x (V x B), (2.11) 



n 

For this system, we have the following magneto-helicity conservation relation: 

— f B-Adx + l\ B-(VxB)dx = 0, (2.12) 

where A such that B = V x A is any vector potential of B. To prove this relation, we 
note, using Green's formula, that 

d 
dt 

Then, taking the scalar product of ( 12. lip with 2 A easily gives the result. 

Remark 2.1 The viscosity term at the right-hand side of ( QE 1 0\) involves the rate of strain 
tensor a (u) = Vm+(Vu) t . In usual gas dynamics, the viscosity term involves the traceless 
rate of strain tensor o q(u) = V«+ (Vu) T — (2/3)(V -u)Id. However, a careful computation 
of the viscosity in the case of isothermal gas dynamics shows that, in this case, the right 
tensor is the full rate- of- strain tensor o~(u) and not its trace-free counterpart o~o(u). 



B ■ A dx = 2 / A- B t dx 
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3 Derivation of the Hall-MHD equations from kinetic 
equations 

In this section, we provide a kinetic formulation of the Hall MHD problem. We start from 
a kinetic equation for the ion distribution function f(x,v,t) of the plasma, where x is 
the position, v the velocity and t the time. This distribution function is a solution of the 
following kinetic equation 

dtf + v- VJ + -{E + v x B) ■ V v f = Q{f) , (3.1) 
m 

where e is the positive ion charge, supposed equal to the absolute value of the elementary 
charge, m is their mass, E(x,t) and B(x,t) are the electric and magnetic fields, and 
Q(f) is the collision operator for electron-ion collisions. We respectively introduce the 
ion density n, mean velocity u and energy W by 

n = [ f dv , nu = j fvdv, W= j fm^—^-dv, 



and the temperature by 



12 



2 f* I u Uj I 

—nT = W n\u\ 2 = / f m dv. 

2 2 1 1 J J 2 

We assume that the electrons are described by their fluid quantities, namely their 
density n e (x, t), their fluid velocity u e (x,t) and their temperature T e (x,t). The use of a 
fluid model for the electrons while the ions are treated kinetically can be justified by the 
small electron to ion mass ratio. A formally rigorous justification of this can be found 
e.g. in [7]. The electron density is supposed equal to the ion density by quasineutrality: 

n e = n . (3.2) 

The electron momentum conservation equation, when the transport term is neglected due 
to their small mass, gives rise to the generalized Ohm's law (see ( 12. 2 j) where all terms in 
factor of m e are set to zero): 

V x {nT e ) + en(E + u e x B) = ennj , (3.3) 

where m e is the electron mass, j is the current density and r) is the resistivity. The current 
density is given by: 

j = en(u — u e ) . (3.4) 

Since the electron mass is neglected and assuming monoatomic gas equation of state, the 
electron energy W e can be expressed in terms of the electron temperature T e by 
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For the same reason, the mass is neglected in the electron energy flux, which reads \nT e u e . 
For simplicity, we consider a model electron-ion collision operator as follows: 

Q(f) = — V„ ■ ((v - u e )f + ^V„/) • (3.5) 
m m 

The first term expresses the relaxation of the ion velocity to the electron one, while the 
second one expresses the relaxation of the ion temperature to the electron one. More 
realistic expressions of the electron-ion collision operator can be found in the literature 
(see e.g. [2]), but this model is chosen for the sake of simplicity of exposition. The rate 
of change of the ion momentum is given by: 

J Q{f) mv dv = —er/nj, 

and is the opposite of the right-hand side of (13. 3p . which is consistent with the total 
momentum conservation of the electron-ion collisions. The rate of change of the ion 
energy is given by 

Q(f) -^r~ dv = -ewj ■ u + 3 [T e - T). 

2 m 

By total energy conservation in electron-ion collisions, the rate of change of the electron 
total energy is the opposite. Then, the electron energy conservation equation reads: 

d t (^riT^j + V x ■ (^nT e u^j = -enu e ■ E + er]nj ■ u + 3 6 ^ (T - T e ). (3.6) 

In this equation, the first term of the right-hand side is the work done by the electrons in 
the Lorentz force, while the last two terms are due to the electron-ion collisions. Taking 
the scalar product of ( 13 . 3 H by u e and subtracting it to (13. 6ft leads to 

(dt + u e ■ V x .) QnT^ + ^-nT e V x ■ u e = V \j\ 2 + 3^^(T - T e ). (3.7) 

The first term of the right-hand side is Joule heating of the electrons, while the second 
term is the electron temperature relaxation to the ion temperature. 
The magnetic field evolves according to the Faraday equation 

d t B + V x x E = , (3.8) 

and the current is linked to the magnetic field by Ampere's law: 

V x xB = voj , (3.9) 

As in the previous section, the displacement current has been neglected. 
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As a summary, the considered kinetic model is as follows: 

dtf + v ■ V x f + -(E + v x B) ■ V v f = ■ {(v - u e )f + —V v f) , (3.10) 

m mm 

V x (nT e ) + en(E + u e x B) = er/nj , (3.11) 

(d t + u e ■ V x ) + ^T e V x ■ u e = 77|i| 2 + 3^-(T - T e ) , (3.12) 

d t B + V x x E = 0. (3.13) 
V x xB = fi j, (3.14) 
j = en(u-u e ), (3.15) 

n = J f dv , nu = J fvdv, -nT — J f m- — - — —dv. (3.16) 

Now, we link this system to Hall-MHD by taking the moments of the ion kinetic 
equation. Integrating ( I3.10P with respect to v after premultiplying it successively by 1, 

I 1 2 

mv or m^j-, we get the following ion mass, momentum and energy balance equations: 

d t n + V x - (nu) = 0, (3.17) 
m (d t (nu) + V x ■ (nu <g> u)) + V x ■ P = en(E + u x B) - er/nj , (3.18) 

2 2 

d t W + V x • (Wu + Pw + q) = en£ • u - er]nj ■ u + 3^-^-(T e - T) , (3.19) 

m 

where P and q are the stress tensor and heat flux vector, given by: 

/Tfl f 
f (v — u) (g) (v — u) dv, q = — / / (v — u)\v — u\ 2 dv . (3.20) 

By combining (13.181) with the generalized Ohm law ( 13. lip , we obtain the total fluid 
momentum balance 

m (^(nu) + V,, • (nu <g> u)J + V, • (P + nT e ld) = jxB. (3.21) 

Using Ampere's equation ( 13.14p . the fact that V x -5 = and the vector identity (V x x 
B) x B = V x • (B ® B) — Va;(|-B| 2 /2), the total fluid momentum balance can be written 
in conservative form: 

d 1 

— (mnu) + V x ■ (mnu ®u B ®B + P tot ) = , (3.22) 

at /j,Q 

where the total pressure tensor P tot is written 

\B\ 2 

P tot = F+(nT e + l -^)ld. (3.23) 

ZflQ 
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Similarly, by adding the energy conservation equations (13.61 ) and A3. 19j) . we get, for the 
total fluid energy W f = W + W e : 

^W f + V x ■ ((W e + V e)u e + Wu + Pit + q) — E ■ j . (3.24) 
The Faraday equation (13.131) implies that 

4 OP-) + —B-(V x xE) = 0. (3.25) 
dt V 2/i / no 

By adding (13.241) and (I3.25P and using Ampere's law (I3.14|) . the total energy (which is 
the sum of the total fluid energy and the magnetic energy) W to t = Wf + satisfies the 
following conservation law: 

JVtat + V, ■ ((W e + p e )u e + Wu + Fu + q + — E x B) = , (3.26) 
at Ho 

The MHD equations are obtained through the closure assumptions that P = nTId, 
q = 0, which can be justified e.g. by a Maxwellian closure, i.e. assuming that / = M njMj T 
with 

n ( m\v — u\ 
M ^ = — -^exp ^ _ 

v m ) 

The Maxwellian closure itself can be justified if ion-ion collisions are strong enough to relax 
the distribution / quickly to M njU ^. However, in many instances, the Maxwellian closure 
is used in spite of not being fully justified. In this case, we obtain the 2-temperature, 
Hall, resistive compressible MHD equations which are as follows: 

d t n + V x ■ (nu) = , (3.27) 

1 \B\ 2 

d t (mnu) + V x ■ (mnu <g> u B ® B + (n(T e + T) + — )Id) = , (3.28) 

Mo 2/i 

/ 3 \B\ 2 \ 5 1 

d t [W+ -nT e + — + V x • (Wu + nTu + -nT e u e + — E x B) = , (3.29) 
V 2 2/i / 2 fio 

d t B + V x x E = , (3.30) 

V x (nT e ) + en(E + u e x B) = erjnj , (3.31) 

(dt + u e ■ V x ) (~nT e ) + \nT e V x ■ u e = v\j\ 2 + 3^^(T - T e ) , (3.32) 

V x xB = fi j , (3.33) 
j = en(u — u e ). (3.34) 

The first four equations are the basic conservation laws of mass, total momentum, total 
energy and magnetic field. Eq. (I3.3ip is the generalized Ohm law and provides the 
expression for E. Eq. (13.321) provides the evolution equation for T e . Finally, Ampere's eq. 
(13.331) defines j and (13.341) defines u e . The fact that u e ^ u gives rise to the Hall term. 
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Now, if T = T e , then, eq. (I3.32p can be removed and one gets the single temperature 
Hall resistive compressible MHD equations. Another simplification is to suppose that the 
common ion and electron temperatures are constant (isothermal assumption). In this 
case, the total energy equation (13.291) is a consequence of the Faraday and momentum 
eqs. (I3.30p . (13.281) and can be removed. Then, we find the model of section [2j 

A kinetic formulation of the model considered in section [2] is also obtained from the 
kinetic model (I3.10p - (l3.16p by supposing that the electron and ion temperatures are the 
same constant T. This kinetic model is written below: 

2 rp 

dtf + v ■ V x f + -{E + v x B) ■ VJ = ■ ((« - u e )f + -V„/) , (3.35) 

m mm 

TV x n + en(E + u e x B) = erjnj , (3.36) 

d t B + V x x E = , (3.37) 

V x xB = fi j , (3.38) 

j = en(u — M e ) . (3.39) 

By eliminating E, j and u e respectively using (I3.36p . (I3.38p . (I3.39p . and after some easy 
algebra, we find the following coupled Fokker-Planck Faraday system for / and B: 



d t f + v-VJ + - 
m 



1 T 

u) x B H (V x x B) x B V x lnn 

«oen e 

,]n V v -((v-u)f + -V v f), (3.40) 



m m 



d t B + W x x [Bxu + (V x . x B) x B ] = -— V,, x (V^ x B) . (3.41) 

For the kinetic eq. f)3.40p . we have the following entropy dissipation identity associated 
to the entropy H(f): 



dt T J( x ,v)eM. 3 xM. 3 f 
H(f) = [ f In fdxdv. 



1 (x,v)eK. 3 xl 

Particle-in-Cell simulations of this model can be found in [6]. Now, taking the first two 
moments of (I3.40p . we get the isothermal resistive Hall MHD model which was written 
in dimensionless form at (j2.9|) - (j2.11j) (taking v = and p(n) = 2Tn). It is interesting 
to note that in order to get a kinetic model for the standard resistive isothermal MHD 
equations, we need to neglect the Hall term j^;(V x x B) x B in the Faraday eq. (13.41 p . 
but we must retain the corresponding term in the kinetic equation (I3.40p . Therefore, a 
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kinetic formulation for standard resistive isothermal MHD is as follows: 



d t f + v-V x f + 



m 



u) x B 



d t B + V x . x (B x u) = -—V x x (V x x B) . 



T 

-(Va, x B) x B V x lnn 



e 2 rjn 
m 



V, • {{v 



T, 



u)f + -V v f) 
m 



(3.42) 
(3.43) 



We can also get a kinetic formulation for the ideal isothermal MHD equations by neglecting 
the resistive term -^'Vx x (V x x B) at the right-hand side of the Faraday eq. (I3.43p . 
The collision term at the right-hand side of (I3.42p has no contribution in the mass and 
momentum conservation equations and can be either kept or removed without modifying 
the corresponding balance equations. 



4 Existence result for the incompressible, viscous, re- 
sistive Hall-MHD 

In the paper, we use the function spaces 

H per (Vt) = {B E (H 1 ^)) 3 | V • B = on Q with periodicity conditions}. 
We use (A, B) = J n A ■ B dx, for any pair A,B<E (L 2 (f2)) 3 . We remark that 

ll V5 ll(L2(Q))9 = ll V x B W(LHn))3 ■ ( 41 ) 

To prove theorem II. 1[ we just focus on the Hall problem itself. The extension to the full 
Hall-MHD problem is explained at the end of the section. So, our goal is now to prove 
the following existence theorem for the Hall problem. We introduce the following weak 
formulations of the Hall MHD problem: 

Find B e L°°(0, T, L 2 (tt)) n L 2 ((0,T), H per {ti)) such that for any A e H per {ty: 

(A, dtB) + (V x A, (V x B) x B) + (V x A, V x B) = 0. (4.2) 

Theorem 4.1 Assume that Bq G (L 2 (f2)) 3 with V • Bo — 0. Then, there exists a global 
weak solution B for the Hall problem \1.7\) , U.8\) . Moreover, we have B e L°°((0, T), L 2 (Q)) 
D L 2 ((0,T), H per {n)) and d t B G L%((0,T), H~ 2 (Q)). 

The proof is based on the construction of an approximate solution by Galerkin's 
method. Uniform a priori bounds on these approximate solutions will allow us to pass to 
the limit thanks to standard compactness arguments. 

We define the Fourier basis (4>k)kez 3 , with (j) k = e 2mk ' x . We denote by \k\ = \k±\ + \k 2 \ + 
\k 3 \. The Fourier basis provides a complete ortho- normal basis of L 2 (Q). We denote by 
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P the projection onto the Fourier basis with index |i| < N. We define an approximate 
solution B N : (0,T) ^ H per (Q) of problem (fTTT|) . (TOj) of the form: 

ieZ 3 , |i|</V 

with the divergence free constraint i ■ Bf(t) = 0, satisfying 

(A N , d t B N {t)) + <V x A N , (V x B N (t)) x + (V x A N , V x = 0, (4.4) 

for all A N of the form 

TV 

i=i 

and with initial condition B N (0) = P n Bq. 

Lemma 4.2 There exists a global in time solution B N (t) of ^4-4\) which is uniformly 
(independently of N) bounded in L°°((0, T), L 2 {VL)) n L 2 ((0, T), H l {Vt)) such that d t B N 
is uniformly bounded in Li ((0, T), H~ 2 (Q)). 

Proof. Throughout of the proof, C denotes a generic constant. We take A = B N in 
P~4"j) and get: 

(B N ,d t B N ) + \\VxB N \\f L2m3 =0. (4.5) 

Thus, we have: 

\\B N (t)\\l Hn) + 2 [ t \\VB N (s)\\l, m9 ds^ |K||| 2(n) . (4.6) 

This shows the uniform bound of B N in L°°((0, T), L 2 (VL)) n L 2 ((0, T), ^(fi)). 

To show the time regularity, we take A G H 2 (Q) 3 , such that V • A = and take 
a jv = p n a 

(A, d t B N ) = (A N , d t B N ) = - (V x A N , (V x B") x B N ) -(Vx/,Vx B") . (4.7) 
We recall that 

(VxB)xB = V-(B®5)-V(^). (4.8) 

Introducing ( 14. 8ft into (14. 7p and noting that, by Green's formula, a gradient and a curl 
are orthogonal in L 2 , we find 

(A, d t B N ) = (VxA ff ,V' (£" ®B N )) + {W x A N , V x B N ) . (4.9) 

To estimate the first term at the right-hand side of (14. 9p . we use Green's formula and 

get: 

(V x A N , V • (B N <g> = - (B N , (B N ■ V)V x A N ) . 

15 



Using Holder's inequality, and remarking that \\A \\H 2 (n) < ||^4||_f/ 2 (n), we get 



? JV||2 



\(B N ,(B N .V)VxA N )\<C\\B } ma)l 
By Gagliardo-Nirenberg's inequality [TO] . 



.4 



N 



\H 2 (n) 



< C\\B 



N\\2 



lL 4 (H)l 



Allien). (4-10) 



l|5 JV |li^)<C||V^||| 2(n) (||^||| 2(n) + l), 
and using the uniform L 2 bound, we can write: 

| (B N , (B N • V)V x A N ) | < C\\VB N \\l 2{Q) \\A\\ H , {n) . 
The second term is simply estimated by using Cauchy-Schwartz inequality: 

(Vxi ff ,Vx B N ) | <C\\VB N \\ LHa) II^H^n) < C\\VB N \\ L2{a) \\A\\ mm 
Collecting these estimates, we obtain: 



(A, d t B N ) | < C 
< C 



\VB 



N 



\VB 



N 



Il 2 (q) 



\L 2 (fl) 



\A\ 



H 2 (fl)- 



Therefore, 



Thus: 



\9tB \\h- 2 (q) < C 



\d t B \\ 3 H -2 {n) < C 



\VB 



3 

JV|| 2 



\L 2 (Q) 



Thanks to (14. 6p . the right-hand side is integrable on (0,T) and we get 

cT 



\d t B N \\l_ HQ) dt < C, 



(4.11) 



which ends the proof of Lemma 14.21 



Thanks to this Lemma, we can proceed to the 

Proof of Theorem 14.11 According to Lemma [4.21 the sequence (B n )ne® is uniformly 
bounded in L 2 ((0, T), H 1 ^)) with (d t B N ) Nm uniformly bounded in Lt((0,T), H~ 2 (tt)). 



Consequently, by virtue of Lions- Aubin Lemma [T7] , (B 



is compact in L 2 ((0, T), L 2 (fi)). 



Therefore, there exists a subsequence (B Nk )km and a function B in £ 2 ((0, T), if 1 (f2)) with 
d t B in Ll((0,T),#- 2 (ft)), such that: 



B Nk -± B weak star in L°°((0, T), L 2 (0)) 

B Nk -± B weakly in L 2 ((0, T),H l (n)) 

B Nk B strongly in L 2 ((0, T), L 2 ({1)) 

k 3*5^ d t B weakly in L^((0,T), H~ 2 (tt)) 



(4.12) 
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We take A e H 3 (tt) 3 , such that V • A = 0, take A Nk = P Nk A in f Oj) and integrate it 
with respect to time. We have: 

(A N \B Nk (t)) - (A Nk ,B Nk (0)) + f (V x A N \ (V x B Nk ) x ds + 

+ / (V x A* fc ,V x B^*) ds = 0. 

Thanks to the Sobolev imbedding in dimension 3, VxAe L°°(f2) 3 and V x A Nk — > 
V x A strongly in L°°(f2) 3 . Now, thanks to the convergences (14.121) . we can take the 
limit Nk — > oo and get that B satisfies 

(A, B(t)) - {A, B(0)) + / (V x A, (V x B) x B) ds + 

Jo 

+ / (V x A, V x B) ds = 0, 

which is a weak solution of the Hall problem. ■ 

Proof of theorem 11.11 We apply the same Galerkin construction (u N , B N ) for the 
coupled system (jl.ip . (11.31) as we did for the Hall system and use the energy identity for 
the Galerkin approximation: 



\u N (t)\\Uu } + \\B N (t)\\ 2 L2{n) + 2 [\\\Vu N (s)\\l 

Jo 



+ \\^B N ( s )\\ 2 (L 2 (n))9) ds 

- \\,, N \\ 2 J- II R^ll 2 

- \\ u o \\L 2 (n) + ll-°o IIl 2 (q)- 



Then the same proof can be reproduced. We just note that the time regularity of u can 
be improved because there is no Hall term in the velocity equation. Therefore, we find 
u t e Lt((0,T),i7^ 1 (fi)) while B t E Ll ((0, T), ff" 2 (fi)). ■ 



5 Axisymmetric flows, KMC waves and Maxwell reg- 
ularization of the non-resistive Hall problem 

In this section, we assume axisymmetric B field. Let x coordinate system, 

where x is the symmetry axis. Axisymmetry about the x-axis means that, given any 
rotation R about this axis, the field B satisfies: B(Rx) = RB(x). We use (x,r,6) 
the cylindrical coordinates of x and (e x , e r , eg) as associated local basis. Then, using a 
representation of axisymmetric, divergence-free fields given in [18], we can write 

B = be e + V x (^e e ), 



17 



where the scalar functions b and ip are functions of (x,r). We note the simple formulas 

V x B = -Cijjee + V x (bee), £ = V^ r) - i 

V x (V x B) = -Lb + V x (je e 

Vx((Vx5)xB) = 




where the Poisson Bracket {a,b} xr = d x ad r b — d x bd r a. 

In these coordinates, the Hall problem (jl.7p . f ll.Sp is written: 

6, + j£,r^J- j^r&J = Cb. 

If initially ?/>o = (which means that the B field is purely swirling), then ip = at all 
times and the equation reduces to the simple viscous Burger's equation 

2 

b t bb x = £b. 

r 

In the inviscid case, the Burger's equation has solutions in the form of propagating shock 
waves. In the context of Hall MHD, these waves are known as the KMC waves, for 
Kingsep, Mokhov and Chukbar [16] (see also (3]). They only exist if the Hall term is 
present. 

Here, we focus on the non-resistive Hall problem and formulate it as a limit of a 
so-called Maxwell regularization. The non-resistive Hall problem is written 

d t B + V x E = 0, 
VxB = j, 
E = j x B. 

In the case of axisymmetric purely swirling I?-fields, this yields the inviscid Burger's 
equation as seen above. We now consider a regularization of this problem by restoring 
the displacement current in the Ampere equation. This yields the problem 

d t B + V xE = 0, (5.1) 
-sd t E + V x B =j, (5.2) 
E=jxB, (5.3) 

where e <C 1 is the Maxwell regularization parameter. We now investigate two formula- 
tions of the Maxwell-regularized problem. 
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1. The (B,j) formulation. This formulation consists in classically eliminating E for j 
and gives 

d t B + V x (j x B) = 0, 
-£<9 t (j xB) + VxB = j. 

When e — > 0, it clearly tends to the non-resistive Hall problem. 

2. The (B, E) formulation. This formulation consists in eliminating j for E. Prom (15. 3p . 
we can write 

B >< & ^ 

3 = ^bT + xb > 

where A is a priori unknown. However, we also have the geometric constraint 

E ■ B = 0, (5.4) 
and A can be viewed as a Lagrange multiplier of this constraint. Indeed, from 

-ed t E + VxB= + XB - 

\B\ 

and taking the scalar product of this equation with B, we get 

-e{d t E) • 5 + (V x B)-B = X\B\ 2 . 
But, differentiating the constraint (15. 4p gives 

(d t E) ■ B = -(d t B) ■ E = (V x E) ■ E. 

Therefore, 

A = ^(-^(V x E) ■ E + (V x B) ■ B). 

Finally, the Maxwell-regularized Hall problem in the (E, B) formulation reads 
d t B + V x E = 0, 

-ed t E + VxB = r^{B x E + (-e (V x E) ■ E + {V x B) ■ B) B). 
\B\ Z 

In this system, the constraint (15 .4p is satisfied as soon as it is satisfied at t = 0. The fact 
that the limit e — > of this problem leads to the Hall problem is no more so obvious. 
Indeed, the limit e — > leads to the following problem: 

d t B + V x E = 0, 

VxB = -—r(B x E + ((V x B) ■ B) B), 
\B\ 

d t (E-B) = 0. 
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The second equation is equivalent to 

B x [(V x B) x B - E] = 0. 
If (E ■ B)\ t=0 = 0, then E ■ B = for all times; Then, we can invert this equation into 

E = (V x B) x B, 

and recover the Hall problem as the formal limit of the Maxwell-regularized system. It is 
an interesting problem to determine if this limit can be made rigorous. 

6 Conclusion 

In this paper, we have derived and analyzed the Hall-MHD model. First, the model has 
been derived from a scaling limit of a two-fluids Euler- Maxwell system, under suitable 
scaling assumptions. Then, a derivation of the Hall-MHD from a kinetic model consisting 
of a Fokker-Planck equation for the ions coupled with fluid electrons has been realized. 
In the analysis section, we have proved the existence of global weak solutions for the 
incompressible viscous resistive Hall-MHD problem. The proof relies strongly on the skew- 
symmetric structure of the Hall term, which does not affect the energy estimate. This work 
shows that maintaining this structure is crucial for the well-posedness of the problem and 
is likely to be crucial as well for the derivation of stable numerical approximations. The 
question of the perfectly conducting boundary condition will be also investigated in future 
work. In the last section, we have reviewed some aspects of axisymmetric, purely swirling 
magnetic fields and proposed a regularization of the Hall problem by reintroducing the 
displacement current in Ampere's equation. Future work will be devoted to the analysis 
of the axisymmetric case and to the Maxwell regularization of the Hall problem, as well 
as to the investigation of the kinetic formulations of the Hall-MHD problem. 
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